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Effect  of  diurnal  convection  on  trapped  thermal  plasma 
in  the  outer  plasmasphere 


1.  Introduction 

The  effects  of  magnetospheric  convection,  and  its  associated  large-scale  electric  field,  on 
the  global  structure  and  morphology  of  the  magnetosphere  are  well  known  [Axford,  1969]. 
The  most  striking  example  in  the  inner  magnetosphere  is  the  formation  of  the  plasmapause, 
whose  location  is  determined  by  the  interplay  between  the  large-scale  electric  field  and  the 
corotation  electric  field  [Nishida,  1966;  Lemaire,  1974].  The  effects  of  convection  on  the 
detailed  structure  of  the  region  inside  the  plasmapause  are  less  well  known. 

Recent  work  on  plasmaspheric  convection  has  utilized  fluid  methods  [Rasmussen  and 
Schurik,  1 990;  Rasmussen  etal.,  1993;  Khazanov  et  al.,  1994],  which  treat  the  entire  distribution 
of  particles  identically.  This  is  reasonable  if  one  assumes  that  the  flux  tubes  are  completely 
filled,  i.e.,  a  diffusive  equilibrium  condition  exists.  However,  convection  (the  diurnal  motion 
of  a  magnetic  flux  tube)  affects  different  regions  of  the  plasma  phase  space  in  fundamentally 
different  ways,  leading  to  a  steady  state  in  the  high-altitude  regions  near  the  plasmapause, 
but  not  hydrodynamic  or  thermal  equilibrium.  Because  this  partition  in  phase  space  is  not 
realizable  in  a  fluid  formulation,  it  is  necessary  to  treat  the  plasma  from  a  kinetic  point  of 
view,  particularly  at  high  altitudes  where  Coulomb  collisions  are  infrequent.  We  also  note 
that  differences  between  measured  densities  and  fluid  model  predictions  have  been  previously 
reported  [Rasmussen  andSchunk,  1990;  Craven  etal.,  1997]. 

In  this  Letter,  we  use  collisionless  kinetic  theory  to  investigate  the  convection  of  mag¬ 
netically  trapped  thermal  particles,  and  to  elucidate  their  role  in  determining  the  density  and 
temperature  morphology  inside  the  plasmapause.  This  trapped  population,  found  in  high- 
altitude  flux  tubes,  is  thermodynamically  isolated  from  the  ionospheric  plasma  on  transport 
time  scales,  and  is  strongly  affected  by  the  convection  electric  fields  which  form  the  plasma- 
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pause.  These  flux  tubes  experience  expansion  and  contraction  due  to  radial  excursions  during 
their  diurnal  trajectory.  This  affects  the  phase-space  density  of  the  trapped  population.  Flux- 
tube  motion  does  not,  however,  modify  the  untrapped  population  because  those  particles  escape 
and  re-enter  the  ionosphere  on  a  time  scale  (~  1  hour)  [Lemaire,  1989]  shorter  than  that  of 
convection  (~  12  hours). 

In  addition,  the  present  approach  has  the  potential  for  investigating  a  wide  variety  of 
instabilities  (due  to  the  nonequilibrium  nature  of  the  distribution  function)  with  a  possibility 
of  explaining  the  observed  asymmetric  wave  activity  in  the  outer  plasmasphere  [. Boardsen 
et  al.,  1995].  The  waves  could  provide  a  mechanism  for  scattering  particles  between  the 
trapped  population  and  the  untrapped  population.  These  velocity-space  diffusion  effects  are 
not  included  in  the  present  analysis. 

2.  Convection  model 

We  divide  the  plasmasphere  into  two  regions:  the  barosphere,  where  Coulomb  collisions 
dominate  and  the  plasma  distribution  function  is  Maxwellian;  and  the  region  above  the  baro¬ 
sphere,  which  is  taken  to  be  collisionless,  and  where  the  convection  of  trapped  particles  is 
included.  The  boundary  between  these  two  regions  is  the  baropause,  approximated  by  an 
infinitesimally  thin  surface.  Beyond  the  plasmapause,  the  trapped  particles  are  absent  because 
they  are  continually  convected  away  toward  the  magnetopause. 

For  this  study,  we  consider  the  low-energy  limit,  ignoring  both  curvature  and  gradient 
drifts:  all  particle  motion  is  field  aligned,  with  both  energy  and  magnetic  moment  conserved. 
In  this  limit,  there  are  four  classes  of  particles  possible  in  a  closed-field-line  region,  and  they 
can  be  identified  by  their  position  in  velocity  space  [Lemaire  and  Scherer,  1970]:  (1)  escaping 
particles  which  have  small  pitch  angles,  enough  energy  to  escape  the  gravitational  trap,  and  can 
enter  the  other  hemisphere,  (2)  incoming  particles  which  escape  from  one  hemisphere  and  enter 
the  other,  (3)  ballistic  particles  which  do  not  have  enough  energy  to  overcome  gravity  and  return 
to  the  ionosphere,  and  (4)  trapped  particles  which  are  reflected  between  two  mirror  points  and 
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never  encounter  the  baropause.  The  first  three  classes  are  in  contact  with  the  ionosphere  and 
together  make  up  the  portion  of  the  distribution  called  the  source  cone.  Because  of  this  thermal 
contact,  the  characteristics  of  the  source-cone  particles,  such  as  temperature  and  density,  are 
determined  by  ionospheric  conditions.  The  fourth  class  (trapped  particles)  forms  a  loss  cone  in 
velocity  space,  and  can  exhibit  a  morphology  significantly  different  from  the  ionospherically 
generated  source  cone,  even  during  magnetically  quiet  times. 

The  velocity  distribution  is  assumed  to  be  Maxwellian  at  the  baropause,  transforms  into  a 
source-cone  distribution  as  it  flows  up  the  field  line,  and  re-enters  the  barosphere.  The  trapped 
particles  in  each  flux  tube,  however,  remain  in  the  flux  tube  indefinitely  (in  the  limit  of  no 
diffusion),  and  convect  under  the  influence  of  the  combined  convection  and  corotation  electric 
fields.  Only  those  trapped  particles  which  never  encounter  the  baropause  during  their  diurnal 
motion  (and  hence  would  be  absorbed)  are  included.  Also,  only  those  trapped  particles  which 
drift  entirely  around  the  Earth  are  included.  The  distribution  of  the  trapped  particles  is  taken 
to  be  Maxwellian  over  the  permitted  region  of  phase  space,  and  their  density,  relative  to  that  of 
the  source  cone,  is  a  free  parameter.  We  do  not,  however,  specify  a  mechanism  for  filling  the 
loss  cone. 

To  determine  the  source-cone  distribution  quantitatively,  the  lower  boundary  condition  of 
a  given  field  line  is  a  normalized,  isotropic  Maxwellian  distribution  fo  (at  the  baropause)  for 
each  species  with  mass  m 

(  m  y/2  \  n>(»io  +  «io) 

■  /o(vo)=nH2dwJ  exp( — 

where  no,  T0,  and  v0  are  the  baropause  density,  temperature,  and  velocity,  respectively.  The 
source-cone  distribution  fs  at  a  higher  altitude  r  on  the  same  field  line  is  determined  by  applying 
Liouville’s  theorem  (conserving  energy  and  magnetic  moment)  and  taking  into  account  the 
accessibility  condition 


/s(v,r)  =  /0(v0) 


v\\  +  2 


U-Uo 


m 


(2) 


3 


where  0  is  the  unit  step  function,  B  and  U  are  the  magnetic  field  and  the  potential  energy  at  r, 
Bo  and  Uo  are  those  quantities  at  the  baropause,  and  vo(v)  is  the  velocity  transformation  (see, 
e.g.,  Eviatar  et  al.  [1964]).  The  potential  energy  U  consists  of  the  gravitational,  electrostatic 
and  centrifugal  potential  energies.  The  ©  function  in  Eq.  (2)  indicates  the  region  of  phase  space 
that  the  source  cone  occupies.  The  density  of  the  source-cone  population  is  determined  by 
integrating  the  distribution  function,  ns  =  /  d?vfs>  and  one  obtains 

ras(r)  =  n0  (e-^  -  VA  e-v,/0i)  ,  (3) 

where  ^  =  (U  —  Uo)/kT0  is  the  dimensionless  potential  energy  and  A  =  1  -  B/B0.  This  result 
is  equivalent  to  that  obtained  by  Eviatar  et  al.  [1964],  The  first  term  in  Eq.  (3)  corresponds  to 
isothermal  diffusive  equilibrium,  while  the  trapped  particles  are  subtracted  by  the  second  term. 

To  determine  the  trapped  distribution  quantitatively,  we  include  the  effects  of  convection 
on  the  region  of  phase  space  left  undetermined  by  Eq.  (2).  We  assume  that  the  distribution  in 
this  region  is  proportional  to  /0,  but  that  accessibility  limits  the  trapped  distribution  ft  to 


/t(v,r)  =  f?/o(v0) 

x©  (v\  — 

/  R 
X0 


(4) 


Bod-B  L 


vn  +  2 
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vn  -f-  2 


U-Ud 


m 


D 
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>. Bd-B 

where  rj  is  a  free  parameter  which  characterizes  the  density  of  the  trapped  population  relative  to 
the  source  cone,  Bd  and  Ud  are  the  magnetic  field  and  potential  energy  at  the  equator  of  the  flux 
tube  at  its  closest  to  Earth  (which  we  take  to  be  at  dawn — this  is  true  for  a  symmetric  electric 
field),  and  Bod  and  Uod  are  the  same  quantities  at  the  baropause  of  the  dawn  flux  tube.  Those 
parameters  with  a  subscript  d  are  controlled  by  the  electric  field  model,  which  determines  the 
drift  path  of  the  flux  tube.  The  parameter  77  will  depend  on  flux-tube  filling  and  loss  processes,  as 
well  as  the  recent  history  of  the  plasmasphere.  The  first  0  function  in  Eq.  (4)  allows  only  those 
particles  that  do  not  encounter  the  baropause  during  their  diurnal  convection.  That  is,  we  keep 
only  those  trapped  particles  whose  turning  point  is  always  above  the  baropause.  The  second 
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©  function  in  Eq.  (4)  is  the  accessibility  condition  from  dawn  for  convective  motion.  That  is, 
we  keep  only  those  particles  which  actually  drift  completely  around  the  Earth.  Equation  (4) 
is  the  distribution  only  in  the  equatorial  plane;  a  similar  expression  holds  at  higher  latitudes. 
Integrating  over  the  distribution,  the  density  of  the  drifting  trapped  particles,  nt,  is 

nt(r)  =  17710  [\[Md  e~'p,Aoi  ~  \[Ad  ,  (5) 

where  A0i  =  1  -  B/B0d,  Ad  =  1  -  B/Bd,  ij>d  =  (U  -  Ud)/kT0,  and  /?  =  Ad/(  1  -  Ad). 
Figure  1  shows  the  regions  of  velocity  space  which  the  two  distinct  populations  occupy,  at  the 
equator  of  a  typical  flux  tube. 

The  total  density  is  the  sum  of  the  two  populations,  n  =  ns  +  nt,  as  given  by  Eqs.  (3) 
and  (5).  In  these  expressions  the  electrostatic  potential  is  undetermined  and  must  be  calculated 
self-consistently  by  applying  the  quasineutrality  condition  ne  =  to;. 

Besides  the  density,  knowledge  of  the  total  distribution  allows  the  calculation  of  the 
effective  temperature.  The  distribution  functions  calculated  above  exhibit  highly  anisotropic 
temperatures,  which  could  lead  to  strong  instabilities.  The  actual  spatial  regions  of  instability 
will  depend  on  the  the  parameter  77,  as  well  as  the  specific  electric  field  model. 

3.  Equatorial  densities 

To  illustrate  the  result  of  treating  the  trapped  particles  in  the  manner  described,  we  now 
calculate  plasma  densities  in  the  equatorial  plane  using  the  plasmapause  formation  mechanism 
of  Nishida  [1966]  (corotation  and  convection  electric  fields  combine  to  form  a  separatrix),  a 
uniform  dawn-dusk  electric  field  [Kavanagh  et  al.,  1968],  a  dipolar  magnetic  field  aligned 
with  the  rotation  axis,  a  baropause  of  constant  height,  density  and  temperature,  with  only  the 
gravitational  and  electrostatic  potential  energies  included  in  U  (no  centrifugal  potential).  These 
simplifications  allow  us  to  focus  on  the  contribution  of  the  trapped  population,  and  they  can 
easily  be  relaxed.  For  illustrative  purposes,  we  choose  two  ion  species,  90%  hydrogen  and  10% 
helium  at  the  baropause,  along  with  ro  =  3000  km,  To  =  3000  K,  and  a  dawn-dusk  electric 
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field  of  0.25  mV/m.  A  detailed  study  of  parameter  space  is  beyond  the  scope  of  this  Letter. 

Figure  2  shows  the  electron  density  in  the  equatorial  plane  for  -q  =  1.75.  Inside  the 
plasmapause  (shown  with  the  characteristic  teardrop  shape  expected  from  the  electric  field 
model),  the  density  enhancement  in  the  dawn  sector  and  the  density  depletion  in  the  dusk  sector 
are  due  primarily  to  the  convection  of  the  trapped  population.  Outside  the  plasmapause  the 
density  is  low  because  only  source-cone  particles  are  present. 

Because  of  the  importance  of  future  helium  imaging  to  our  understanding  of  the  plasma- 
sphere,  Fig.  3  shows  the  ratio  of  helium  ion  density  to  hydrogen  ion  density  in  the  equatorial 
plane.  This  ratio  varies  rapidly  for  r  £  2 RE,  but  varies  slowly  for  r  ^  2Re>  in  reasonable 
agreement  with  statistical  studies  [Craven  et  al.,  1997],  This  behavior  is  due  to  the  mass 
ratio  of  the  two  species  and  the  interplay  between  the  gravitational  and  electrostatic  poten¬ 
tials.  Figure  3  also  shows  that  this  ratio  is  approximately  constant  across  the  plasmapause,  an 
observation  noted  previously  [Horwitz  et  al,  1986]. 

In  order  to  assess  the  capability  of  imaging  to  distinguish  between  models,  Fig.  4  shows 
the  ratio  of  the  helium  ion  density  given  by  the  present  model  to  what  the  helium  ion  density 
would  be  if  diffusive  equilibrium  conditions  existed.  The  solid  line  indicates  the  surface  where 
the  ratio  is  unity.  Dawnward  of  this  surface,  the  loss  cone  is  overfilled  relative  to  the  source 
cone;  duskward  of  this  surface  (and  outside  the  plasmapause),  the  loss  cone  is  underfilled.  The 
strong  dawn-dusk  asymmetry  of  this  ratio  affirms  the  distinguishability  of  the  present  model. 
To  illustrate  the  effect  of  varying  77,  the  dashed  and  dotted  lines  in  Fig.  4  indicate  the  positions 
of  the  unity  surfaces  for  the  cases  where  77  =  1.5  and  77  =  2.0,  respectively.  (The  densities  for 
these  two  cases  are  not  shown,  but  they  exhibit  a  strong  asymmetry  as  well.) 

A  strength  of  the  present  method  is  that  it  is  not  limited  to  a  particular  electric  field 
specification  or  plasmapause  formation  mechanism;  other  electric  field  models  [e.g.,  Volland, 
1973]  or  plasmapause  models  [e.g.,  Lemaire,  1974]  can  be  used.  In  addition,  it  is  possible  to 
incorporate  the  effects  of  the  penetration  of  the  polar  convection  electric  field  to  low  latitudes 
[Spiro  etal,  1988;  Doe  etal.,  1992]. 
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4.  Conclusion 


We  have  shown  that  trapped  particles  drifting  along  diurnal  convection  equipotentials  and 
adjusting  their  density  self-consistently  leads  to  a  morphology  of  the  outer  plasmasphere  which 
is  significantly  different  from  that  obtained  using  the  assumption  of  diffusive  equilibrium.  The 
density  exhibits  a  strong  asymmetry  between  the  dawn  and  dusk  sectors  compared  with  a 
diffusive  equilibrium  assumption.  The  partition  in  phase  space  between  trapped  particles  and 
untrapped  particles,  and  the  different  effects  of  convection  on  the  different  classes  of  particles, 
requires  a  kinetic  theory. 

Future  spacecraft  missions  will  have  the  capability  to  globally  image  the  inner  mag¬ 
netosphere  (e.g.,  through  scattering  of  the  solar  30.4  nm  line  by  He+),  and  provide  signif¬ 
icant  progress  in  our  understanding  of  the  physical  processes  in  the  plasmasphere  [Meier, 
1991;  Williams  et  al.,  1992],  The  present  model  has  a  predictive  capability  [Reynolds  et  ah, 
1997]  which  is  needed  to  extract  quantitative  information  from  the  observational  data. 

Acknowledgments.  We  would  like  to  acknowledge  fruitful  discussions  with  R.  R.  Meier,  J.  M. 
Rcone,  and  A.  C.  Nicholas.  This  work  is  sponsored  by  the  Office  of  Naval  Research  and  NASA. 
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Figure  1.  Velocity  space  at  the  equator  of  a  typical  flux  tube.  Region  I  is  the  source  cone,  and 
region  HI  is  the  trapped  population.  Particles  in  region  II  encounter  the  baropause,  and  region 
IV  is  not  accessible  from  dawn  due  to  convection.  The  three  boundaries  between  the  regions 
are  determined  quantitatively  by  the  three  ©  functions  in  Eqs.  (2)  and  (4). 
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Figure  2.  Electron  density  ne  in  the  equatorial  plane,  scaled  to  the  electron  density  at  the 
baropause,  neo,  using  the  GSM  coordinate  system.  The  parameters  are  tj  =  1.75,  r0  =  3000 
km,  T0  =  3000  K,  90%  hydrogen  and  10%  helium  at  the  baropause,  and  a  dawn-dusk  electric 
field  of  0.25  mV/m. 
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Figure  3.  Ratio  of  helium  ion  density  nHe  to  hydrogen  ion  density  in  the  equatorial  plane, 
for  the  same  parameters  as  Fig.  2. 
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Figure  4.  Ratio  of  helium  ion  density  nHe  to  diffusive  equilibrium  helium  ion  density  nHe,de 
in  the  equatorial  plane,  for  the  same  parameters  as  Fig.  2.  The  solid  line  indicates  the  surface 
on  which  this  ratio  is  unity.  Near  the  plasmapause,  this  ratio  is  ~  1.5  at  dusk  and  is  ~  0.25  at 
dawn,  a  factor  of  6  difference.  The  dashed  and  dotted  lines  indicate  the  positions  of  the  unity 
surfaces  for  the  cases  where  7/  =  1.5  and  r)  —  2.0,  respectively.  (The  densities  for  these  two 
cases  are  not  shown.) 
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